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Saturation models - colour glass condensate and string percolation - impose a 
strict relation between the average transverse momentum, (-Pt), and the rapidity 
particle densities, dn/dy. By combining this relation with an appropriate evolution 
equation for dn/dy, and imposing energy-momentum conservation, we obtain a fair 
description of data, for generic AB collisions (hadron-hadron, hadron-nucleus and 
nucleus-nucleus) at all rapidities and (high) energies. Predictions are given for the 
LHC. 

1 Introduction 

Over the last few years the relationship between the phenomenological suc- 
cessful String Percolation Model (SPM) [1] and the theoretically sound Colour 
Glass Condensate (CGC) approach [2-12] to high energy hadronic processes 
has been clarified [13-18]. Such a relationship provides not only much needed 
theoretical justification for the assumptions underlying the SPM, but also a 
simplified framework in which to perform phenomenological studies account- 
ing for the essential features of the CGC. 

In [16] it was argued that complementing the SPM by a simple evolution 
equation — which describes the generation of lower rapidity strings from 
higher rapidity ones — yields an effective description of the Glasma [19-24], 



i.e. of the highly coherent state formed immediately after the collision of two 
CGC sheets. Thus, the SPM incarnation of [16] can be used to address the 
early post-collision dynamics in a heavy ion collision. A series of predictions 
[17, 18,25], in which RHIC data was used to fix the model parameters, for 
observables at the forthcoming LHC Pb-Pb and p-Pb programmes have been 
obtained in the framework of this model. Importantly, the results obtained 
for Forward-Backward correlations [18] were found to be in agreement with 
those obtained in the CGC/Glasma [15]. 

In the present work we address the effects due to the inclusion of average 
transverse momentum in our model. From the observation that saturation 
models — CGC and SPM — impose a strict relation between average trans- 
verse momentum and particle density distributions in rapidity, we obtain, 
imposing energy momentum conservation [17], both particle density distri- 
butions and average transverse momentum distributions in rapidity. 

The paper is organized as follows. In Sec. [2] we describe our model. Our 
results, including both a comparison with available RHIC data and predic- 
tions for the LHC, are presented in Sec. O Sec. H] contains our conclusions. 

2 Model setup 

The SPM [1] deals, as its fundamental degrees of freedom, with the longi- 
tudinally extended strings formed in the collision. The relevant parameter 
to address percolation phenomena is the density ^0 which, in the impact 
parameter plane, is defined as 

?-(^)'iv., (1) 

where Ng is the number of strings, or discs in the impact parameter plane, 
r is the radius of each disc, and R is the effective radius of the interacting 
overlapping area. For an AB collision, with A < B, the area of overlap is 
determined by A or, more precisely, by the number of participating nucleons 
from A, i.e. by (1 < ^ A). In the following we shall make the natural 
and simplifying assumption that 

R - RpN'f , (2) 

"'^We will denote the density by ^, instead of the traditional 77 in order to avoid confusion 
with the pseudo-rapidity variable. 
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where Rp is the nucleoli radius. 

The average number of colhsions A'^coi can be estimated given the numbers 
Na and Nb of participating nucleons from each nucleus. As the overlap area 

2/3 

is given by A^^ , the Nb participating nucleons from B are the ones crossing 

1/3 

that area, such that each of them interacts with iV^ nucleons from A. We 
thus have 

Ar,„i ^ n'/Nb . (3) 

When A = B, reduces to the well known Glauber result and for A^^ = 1, 
A'coi = Nb as one would also expect. 

Two important observations are essential for our argument: 

(i) There is a definite relation between the averaged transverse momentum 
{Pt) and the rapidity particle distribution dn/dy [26,27]. The existence of 
such a relation is related to the presence of a single scale in the problem, be 
it the saturation scale in the case of the CGC, or, in the case of Schwinger 
strings, that the field and the charge are related to each other by Gauss' 
theorem. 

(ii) The appearance of low rapidity objects (gluons or strings) arising 
from fast coloured objects (glasma fields or valence strings) can be described 
by an evolution equation [13, 16, 19,20]. 



2.1 The relation between {px) and dn/dy 

The occurrence of percolating behaviour depends on whether strings do or 
do not overlap. If there is no overlap of strings, percolation does not occur 
and one expects the average transverse momentum (pf.) to be given simply 
by the average transverse momentum p\ associated with production from a 
single string, 

= pI and ^ = NsUi , (4) 

where ni is the particle density associated with a single string and A^^ is the 
number of strings. 

On the other hand, when strings overlap, percolation will occur and colour 
vector summation will result both in a suppression of the overall colour charge 
and in an enhancement of the effective string tension, i.e. 

^ = FiONsU, , (5) 
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which can be rewritten, using ([T]) and ([2]), as 



2 



1 * F(f)d^)%, (6) 



jVj/^ dy \ r 

For the average transverse momentum we get 

In both ([6]) and ([7]), is the colour suppression factor given by [28] 




(8) 

The only varying quantity in the right hand side of both and ([7j) is 
the transverse density ^. Thus, a direct relation between I/N"^^ ■ dn/dy and 
(py) necessarily exists [26,27,29,30]. In fact, as F{^) is a decreasing function 
of ^ and F{^)^ is an increasing function of ^, there is a monotomic functional 
dependence between dn/dy and (pj-)- This result is essential in our work. 

Eq. (171) for the average transverse momentum (p^) does not account for 
any kinematical effects. As ^ — 0, F{^) 1 and (I7j) yields (p^) p\ 
(the average transverse momentum of a single string), whereas one would 
expect that in this limit and for a particle of mass m and with rapidity ym = 
In ^Js/m, where ^/s is the centre of mass energy, the averaged transverse 
momentum at the end of phase space, should be zero. Thus, eq. (171) 

needs to be kinematically corrected to the form 

^^^^^W) [i-{soisY] + {s,/syco.h\ ' 

where t] is the pseudo-rapidity. The parameters sq and a in ([ID are fixed 
by the standard relation between rapidity and pseudo-rapidity, taken in the 
limit y ^ y-m and ?7 — cx), to 



such that 



So = 4pf, a=l, (10) 



F(0[l-(4p?/s)] + (4pf/s)cosh2r7" ^ ' 
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The ratio Sq/s is, for high enough energies, very small. Thus, the correction 
dni) to ([7]) is only important for very large values of rj. 

In the high density region ^ 3> 1 corresponding, as we shall see, to mid- 
rapidity, the number of produced strings is proportional to the number of 
collisions, i.e. iV, ~ N^^i [31,32]. Then, from (P, ([2]) and ([3]), it follows that 

e-iV-^/'iVs, (12) 

and that, from ([5]) (or ([6])), 

^ ~ NTe^' - nTn]/' . (13) 

In the low density region ^ <^ 1 corresponding to the forward (backward) 
rapidity region, one has A^^ ~ Na {Ng ~ A''^) [31,32]. Thus, we have 

^ ~ N^J^ (Forward) , (14) 

e ~ N^^^^Nb (Backward) , (15) 

and 

— ^ Na (Forward) , (16) 

ay 

dfi 

— Nb (Backward) . (17) 
dy 

For symmetric collisions Na = Nb, we obtain, as expected, dn/dy ~ Na, 
asymptotically, for all rapidities. However, in the asymmetric case Na Nb, 
flTB]) and flTTl) do not accurately account for the relevant physics. In fact, what 
is implicit in (fT6l) and (fT7|) is that nucleons interact in pairs and that each 
nucleon only interacts once, so that particle production is proportional to 
that from pp with a proportionality factor Na (Nb)- When A^^ < A''^ at 
least some of the nucleons from A reinteract with nucleons from B being 
decelerated relatively to the result of (ITB]) . Correspondingly, in the backward 
direction the nucleons from B are less decelerated than the result from f|T7|) . 
In conclusion, dn/dy is overestimated in f|T6l) and underestimated in f|T7|) . We 
account for this by introducing corrections factors (Na/Nb)'^ and (Nb/Na)'^ 
in ([TB]) and ffT7|) respectively. Thus, we have 

'^'^ Na(^Y (Forward), (18) 



dy \Nb 

dn fN 

Nb{ — 

dy ""KNa 



^r^Ns ] (Backward) . (19) 



5 



Clearly, fll6l) and fll7p are recovered in the limit = Nb- The phe- 
nomenological parameter 7 is to be adjusted. The effect encoded in ( ITSil and 
f fTOj) has been observed [33, 34] in forward production of vr^ in pp collisions 
with 7 ~ 0.11. 



2.2 The evolution equation for dn/dy 

We will use the framework developed in [16] to describe the generation of low 
rapidity strings from a high forward rapidity (valence) string. In a nutshell, 
this generation is given by the logistic equation for population dynamics 



a(-A) s'\ PY, 

where p = p{A,Y) = dn/dy is the particle density, Y = In^y/s/mi,) is the 
beam rapidity and 

A^T^-y, (21) 

with 7] the pseudo-rapidity. The variable —A plays the role of evolution time. 
The parameter 6 controls the low density evolution of p and must therefore 
depend on intrinsic parameters of the theory, py is the asymptotic saturation 
density. 

The solution of (l20l) must be such that 

dp 



both at p = and p = py, and 

d^p 



(22) 



d{-Ay 



= 0, (23) 

Ao 

where the scale Aq defines the separation between the region A > Aq, of low 
density and positive curvuture, and the region A < Aq of high density and 
negative curvuture, and is such that po = p(^0!^) = Py/2. The integration 
of ([20D yields 

p(A,y)= J^l , (24) 
+ 1 

which is nothing but a generalization of the Fermi distribution, known to 
approximately fit RHIC data for pseudo-rapidity distributions in central 
nucleus-nucleus collisions [35,36]. In [36] it was argued that 

py = e^^ , (25) 
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and 

Ao = -«r, (26) 

where a and A are positive constants and a, A < 1. The distribution (12^ is 
then characterized by three constants: a, 5 and A. 

In a colhsion one has not only forward emission, but also backward 
emission. For symmetrical situations, (1241) describes both forward emission 
(?7 > 0) and backward emission {j] < 0). An analogous line of reasoning 
underlies related approaches in the context of the CGC [37]. 

In the most general situation, with Na participant nucleons moving in the 
forward direction and Nb participants moving in the backward direction, the 
fulfilment of ([13]), ([HD and ([HD, leads to 

and 

dn ^ e _ 

Note that in general the forward (F) and backward (B) distributions meet 
at r^, = 5(i + 7)ln^<0. 

If Na = Nb, f l27|) and fl28|) become mirror distributions with, for the 
forward region. 



dn Na e^^ 



dy e"^^ + 1 ■ ^^^^ 
In order to allow for direct comparison with experimental data, the mul- 
tiplicity distributions — ( l29l) for symmetric collisions and ( 1271) and ( l28l) for 
asymmetric ones — must be rewritten as distributions in pseudo-rapidity t]. 
Performing the standard transformation the multiplicity distribution in the 
symmetric case is given by 



1 dn _ e^^ 
For an asymmetric collision we obtain 



— = J — . (30) 



' J ''^-a-.). (F), (31) 



{NaNbY'^ dr, ^Nb/NaY/'^^' e"^^ + 1 
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and 

1 dn e^^ 
{NaNsY/^ Tri = ^ (N^/Nsy/^^^e^^'-'^^l ' ^^^^ 

In both cases the transformation Jacobian is given by 

J ^dy_^ coshr] ^ ^ rri^ +pI 

dv + sinh^ r/ ' v\ 

2.3 Energy- momentum conservation 

As we are deahng with dn/dy distributions for all rapidities, or pseudo- 
rapidities, and as we have some control on the px distribution, via 
it is natural to impose, in some simplified scheme [17], energy-momentum 
conservation. 

We write for the overall energy E conservation 

E= {niT) coshy—dy = (Na + Nb)^ , (34) 

and for the overall longitudinal momentum 

P= {rriT) sinhy^dy = (N^ - N^)^ , (35) 

^ —ym ^ 

where the transverse effective mass (mr) is given by 

{mr) = {pD + m2 , (36) 

with m an effective mass (introduced to avoid the sum over all species). 

In our comparisons with data we impose (IMll and (1351) . Further, in or- 
der to go from the experimentally measured dn/drj (charged) to the dn/dr) 
(total), we use a multiplying factor, consistent with dominance of vr"*", vr", vr" 
production, which is of the order of 3/2. It is, however, possible to construct, 
from (IMl) and fl55|) . a ratio which is independent of multiplying factors 

M.i^^. (37) 
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3 Results: comparison with experimental data 



The model parameters a, 6 and A in the multiphcity distributions for both 
the symmetric case (l30ll and for the asymmetric one, (13T]) and (132|) . along 
with the parameter 7 (only relevant for asymmetric collisions) were fixed as 
to provide the best model joint description of RHIC data (^i = 200 GeV) for 
Au-Au central collisions and d-Au collisions in several centrality classes [38], 
and to further satisfy energy-momentum conservation as given by flMl) and 
fl35|) . The values obtained from this procedure are a = 0.347, 6 = 0.851, 
A = 0.270 and 7 = 0.21. 

In Fig. [1] we show the results obtained from ( l30ll and the experimental 
data for RHIC Au-Au central collisions (^i = 200 GeV) together with our 
prediction for the LHC (Pb-Pb central collisions at ^/s = 5.5 TeV. Here, k 
in the Jacobian (!33|l was set to fixed value k = 1.37. 




Figure 1: l/NAdn/drj in central collisions for Au-Au at ^/s = 200 GeV (full 
line), data from [38], and prediction for LHC Pb-Pb the at y/s = 5.5 TeV 
(dashed line). 



Fig. [2] we plot the multiplicity distributions for several centrality classes 
in d-Au collisions at ^/s = 200 GeV together with RHIC data [38]. The 
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model curves were obtained from ([21]) and (152]) . Here, k (15^ was not set to 
a fixed value, but rather as a function of the average transverse momentum 
Pt ~^ {Pt) given by f|TT]) with m —>■ (m) = 0.23 GeV. The curves fulfil the 
energy-momentum conservation relations flM]) and fl35]) . 




Figure 2: dn/dt] at different centralities from [38] and model curves computed 
from (I3T]) and (132]) . Our curves satisfy energy-momentum conservation in the 
sense of ([34]) and ([35]). 

The average transverse momentum (px) is given, in our model, by fllip . 
The single string average transverse momentum is determined from the curves 
in Fig. [21 whereas the ratio between the string and proton radius r/Rp was 
set to the reasonable value of 1/4. Fig. [3] shows the average transverse mo- 
mentum distribution in pseudo-rapidity for both RHIC and LHC symmetric 
collisions. 

A more detailed comparison between the curves obtained from our model 
and available experimental data both shown in Fig. [2] is given in Table [TJ 
The number of participant nucleons A^^^ and N^^ for each centrality class 
extracted from data [38] is compared with the values used in fl3T]) and f[32]) 
to produce our model curves. The number of collisions Ncoi extracted from 
data is compared with that obtained from ([3]). Further, The experimentally 
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Figure 3: Transverse momentum distributions in pseudo-rapidity (fTT!) for 
central Au-Au collisions at ^/s = 200 GeV and Pb-Pb at ^/s = 5.5 TeV 
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measured total number of particles with |?7| < 5.4 is compared with our 
model prediction computed from the integral of the curve over the same 
pseudo-rapidity interval. 





Nd 




Nau 




A^col 




I? <5.4 


centrality 


(exp) 


(mod) 


(exp) 


(mod) 


(exp) 


(mod) 


(exp) 


(mod) 


0-20% 


2.0 ±0.1 


1.9 


13.5 ± 1.0 


14.5 


14.7 ±0.9 


18.0 


157 ±10 


143 


20-40% 


1.9 ±0.1 


1.8 


8.9 ±0.7 


8.2 


9.8 ±0.7 


10.0 


109 ±7 


106 


40-60% 


1.7±0.2 


1.5 


5.4 ±0.6 


4.8 


5.9 ±0.6 


5.5 


74 ±5 


74 


60-80% 


1.4±0.2 


1.2 


2.9 ±0.5 


2.4 


3.1 ±0.6 


2.6 


46 ±3 


47 


80-100% 


1.1 ±0.2 


1 


1.6 ±0.4 


1.2 


1.7±0.5 


1.2 


28 ±3 


30 



Table 1: Comparison of d-Au {^/s = 200 GeV) data [38] for different central- 
ities with model. A^^^, Nau (mod) are those used in fl3T]) and fl32|) to compute 
the curves shown in Fig. [2l Ncoi (mod) computed with ([3]) using A^^, Nau 
(mod). A'|'5||<5 4 (mod) computed as the area under the model curves in Fig. 

12 

Table [2] shows the explicit fulfilment of energy-momentum conservation 
(!37|) . The LHS is computed from and (l35l) using the model predicted 
average transverse momentum {p'^) ([9]) with as before m (pt) = 0.23 GeV. 
The RHS is computed from the values of Nd, Nau (mod). 





eq. 


(EZD 


centralities 


LHS 


RHS 


0-20% 


0.76 


0.77 


20-40% 


0.63 


0.64 


40-60% 


0.52 


0.52 


60-80% 


0.33 


0.33 


80-100% 


0.09 


0.09 



Table 2: Test of energy-momentum conservation (13 7p 

Finally, predictions for both the multiplicity distribution and average 
transverse momentum for LHC p-Pb collisions {^/s = 8.8 TeV) are shown in 
Fig. H 
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Figure 4: Predictions for multiplicity (left) and average transverse momen- 
tum (right) distributions as obtained from, respectively, (131]) and (132|) . and 

m 



4 Conclusions 

In this paper we were able to present a global description of particle densities 
and average transverse momentum (px) for both symmetric and asymmetric 
collisions for all rapidities, centralities and (high) energies. Although our 
work was carried out explicitly within percolation theory, we believe that 
equivalent results should be obtainable in the CGC framework. There, how- 
ever, the explicit dependence on the number of participants A^^ and Nb, on 
the number of collisions Ncoi, and control of energy- momentum conservation 
have not been worked out in detail. 

Some of our assumptions still lack full theoretical justification. Namely: 

• The relation for 1/NAdn/dy, f l29|) was obtained by taking the limit 

^ oo in the dense region, at mid rapidity, and the limit ,^ — > in 
the dilute region, at large rapidity. Corrections should be expected at 
finite Y and rj. The same apllies to the asymmetric formulae (1271) and 

m- 

• The quantities {Na/NbV and {Nb/NaV, in ([27]) and ([28]), representing 
rescattering corrections necessary if Na 7^ Nb-, are purely phenomeno- 
logical without a theoretical ground (Glauber calculus for instance). 
We do not even know why 7 (forward) = 7 (backward) as suggested 
by the comparison with data. 
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• Our formula for (pl^), ([2]), describing correctly the kinematical (p|^)?7 
correlation in the r] — > ±00 limits, and making finite and meaningfuU 
the integrals flMl) and fl35|) . is presumably not unique. We noticed that 
almost 100% of the integrals flM|) and fl35l) are contained in an interval 
[—1], +7]] equal to [—Y, +Y] where Y is the beam rapidity and thus the 
use of one or another asymptotic formula is not so critical. 

These issues clearly warrant future work. 
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